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Interelectron interactions and correlations in quantum 
dots can lead to spontaneous symmetry breaking of the 
self-consistent mean field resulting in formation of Wigner 
molecules. With the use of spin-and-space unrestricted 
Hartree-Fock (sS-UHF) calculations, such symmetry break- 
ing is discussed for field-free conditions, as well as under the 
influence of an external magnetic field. Using as paradigms 
impurity-doped (as well as the limiting case of clean) two- 
electron quantum dots (which are analogs to helium-like 
atoms), it is shown that the interplay between the interelec- 
tron repulsion and the electronic zero-point kinetic energy 
leads, for a broad range of impurity parameters, to forma- 
tion of a singlet ground-state electron molecule, reminiscent 
of the molecular picture of doubly-excited helium. Compar- 
ative analysis of the conditional probability distributions for 
the sS-UHF and the exact solutions for the ground state of 
two interacting electrons in a clean parabolic quantum dot 
reveals that both of them describe formation of an electron 
molecule with similar characteristics. The self-consistent field 
associated with the triplet excited state of the two-electron 
quantum dot (clean as well as impurity-doped) exhibits sym- 
metry breaking of the Jahn- Teller type, similar to that under- 
lying formation of nonspherical open-shell nuclei and metal 
clusters. Furthermore, impurity and/or magnetic-field effects 
can be used to achieve controlled manipulation of the forma- 
tion and pinning of the discrete orientations of the Wigner 
molecules. Impurity effects are futher illustrated for the case 
of a quantum dot with more than two electrons. 
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I. INTRODUCTION 

Two-dimensional (2D) quantum dots (QD's) created at 
semiconductor interfaces with refined control of their size, 
shape and number of electrons are often referred! - !] to as 
"artificial atoms" . This analogy suggests that the physics 
of electrons in such man-made nanostructures is closebz 
related to that underlying the traditional descriptionel 
of natural atoms (pertaining particularly to electronic 
shells,-and the Aufbau principle), where the electrons are 
takenO to be moving in a spherically averaged effective 
central mean field (CMF). However, using as paradigms 
impurity-doped (as well as the limiting case of clean) two- 
electron-QD (2e QD) analogs to He-like atoms, we show 
that the interplay between the interelectron repulsion 



(Q) and the electronic zero-point kinetic energy (K) may 
lead, for a broad range of impurity parameters, to spon- 
taneous symmetry-breaking (SB) of the self-consistent 
mean field, resulting at zero magnetic field (B = 0) in for- 
mation of a singlet ground-state electron molecule. Such 
SB is beyond the CMF picture and, while negligible in 
the ground state of the He atom (whose study was central 
to the development of the quantum theory of matter due 
to the failure of the Bohr- type modelsa), it is similar in 
nature to the SB foundQ in the 1970's in doubly excited 
He, where formation of an e-He 2+ -e "triatomic" molecule 
has been invoked. Furthermore, we show that symmetry 
breaking at B = of the self-consistent field associated 
with the triplet excited state of the 2e QD originates 
from a Jahn- Teller distortion of the CMF, similar to tl 
underlying formatkwa-of nonspherical open-shell nuclei! 
and metal clusters.BEJJ 

Along with a unification of concepts pertaining to 
spontaneous SB in a variety of finite fermionic systems 
(from nuclei, metal clusters and natural atoms, to 2D 
QD's), we demonstrate the ability to control the orien- 
tation and to manipulate (i.e., to enhance, but also to 
counteract and even to void) the formation of the electron 
molecules in 2e QD's via impurity and/or magnetic-field 
effects. The ability to control the orientation of the elec- 
tron molecule may in principle open new possibilities for 
designing "on-off" (i.e., switching between two discrete 
states) devices, which eventually may be employed in ap- 
plications of QD's as nanoscale logic gatestJ (the effect of 
impurities on the structure of multi-electron molecules in 
QD's with more than two electrons is further illustrated 
in Appendix A). 

That electrons in extended media may undergo crys- 
tallization at low densities, when Q dominates over K, 
has been predictecO by Wigner in 1934. Such Wigner 
crystallization (WC) in. clean QD's results in formation 
of electron moleculesE^TliS [also referred to as Wigner 
molecules (WM's)], which are associated with sponta- 
neous SB, where the symmetry of the ground state, cal- 
culated at the mean-field (self-consistent-potential) level, 
is found to be lower thatuthat of the exact hamiltoaian 
describing the systepiOilj In clean QD's, formationEj of 
WM's is controlledli3£y by the parameter Rw — Q/K. 
For a parabolic confinement (with frequency loq) at B = 
0, it is customary to take Q = e 2 / kIq and K = Tiujq, 
where k is the dielectric constant and lo = {Ti/m*L)o) 1 / 2 
is the spatial extent of the lowest state's wave function, 
of an electron with an effective mass m*; WM's occurEEl 
for Rw > 1, corresponding to much higher electron 
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densitieal3'ta_Jhan those predicted for WC in an infinite 
2D medium^ 

The many-body hamiltonian for a QD with N e elec- 
trons can be expressed as a sum of a single-particle part 
and the two-particle interelectron Coulomb repulsion, 

JV e JV e Ne 2 

I— 1 2 — 1 J>2 J 

The contributions to the single-particle part are written 
as 

H(i) =H (i)+H B (i) + Vi(i) , (2) 

and they contain a term describing the motion of an 
electron in a 2D parabolic confinement, i.e., Hq(i) = 
pf/2m* + m*ujQ(xf + yf)/2, where cj is the frequency of 
the 2D isotropic harmonic confining potential. Magnetic- 
field effects are included in Hb(i) = [(Pi — eAj/c) 2 — 
pf]/2m*+g*^sB-Si/?i, where the vector potential, Aj = 
B(~y i /2,Xi/2,0), is taken in the symmetric gauge, and 
the last term is the Zeeman interaction with an effective 
factor g* , Si is the electron spin, and /j,b the Bohr mag- 
neton. To include the effect of (Coulombic) impurities, 
we added to H(i) the term Vi{i) = {e/n) J2i Qi/\ r i — T^-i\ 
(i = 1, N e ), where Qi is the charge of the Zth impurity 
located at TZi = {xi,yi,di); such impurities which in gen- 
eral may be situated out of the 2D plane of the QD (that 
is di ^ 0) may correspond to implanted atoms (donors or 
acceptors) or represent a fabricated, controllable voltage- 
gate. 

A clean QD (that is with Qi — in the above hamil- 
tonian) may be regarded as a realization of the Thomson 
atom, as are jellium models of metal clusters,E£l where 
the positive charge is uniformly distributed; in the 2e 
Thomson QD (TQD), the confinement to the 2D plane 
is modeled by a parabolic potential. On the other hand, a 
QD analog of the Rutherford model of the atom (RQD), 
where the positive charges are grouped in a single nu- 
cleus, can be achieved through the introduction of a cen- 
tral attractive impurity (in addition to the harmonic con- 
finement). 

The presence of many free parameters in Eqs. (|l|) and 
results in a plethora of possible case studies generated 
by varying the material dependent parameters (k and 
to*), the harmonic confinement loq and the magnetic-field 
B, as well as the number of impurities Z, their charges Qi, 
and positions IZi. In this paper, we have chosen to discuss 
here the following three representative classes of cases at 
a specific value of Rw — 2-39, which is sufficientlly high 
so that the electrons formll3 a WM in the case of a clean 
QD, thus allowing for systematic investigations of the 
effects of impurities on the formation, orientational pin- 
ning, and structural distortions of the electron molecules 
(for the otherparameter values used throughout this pa- 
per, see Ref. |23|). The three representative classes, which 
we discuss, are: (i) A 2e TQD at B — in the presence of 
two off-centered impurities situated on both sides of the 



dot (section II); (ii) A 2e RQD (with a single attractive 
central impurity) in the presence of an applied magnetic 
field (for both the cases of a weak and a strong field, sec- 
tion III); and (iii) A QD with eight electrons at B = 
and a single impurity with varying strength, polarity, and 
location (Appendix A). 

For the case of a clean 2e QD, with the hamiltonian 
given by Eqs. (0) and (§) with N e = 9_apd Vi = 0, the 
exact solution can be found quite easilycilEj owing to the 
separability of the Schrodingcr equation in the center-of- 
mass and relative coordinates; in the presence of impu-. 
rities, separability is lost and finding an exact solutionis 
becomes significantly more complicated (even for two 
electrons in the presence of a central impurity). There- 
fore, and also in order to expound the physical princi- 
ples underlying spontaneous symmetry breaking in QD's, 
we will mostly use in the following (sections II and III 
and Appendix A) the self-consistent spin-and-space un- 
restricted Hartree-Fock (sS-UHF) methcdEfLwhich, un- 
like the restricted HF (RHF) technique ,E2rE3 allows for 
the formation of broken-symmetry states (of lower en- 
ergy than those obtained via the RHF). This sS-UHF, 
which we introduced for studies of QD's in Ref. [l6L em- 
ploys N e (where N e is the number of electrons) orbital- 
dependent, effective (mean- field) potentials and it differs 
from the usuaEj RHF in two ways: (i) it relaxes the 
double-occupancy requirement, namely, it employs dif- 
ferent spatial orbitals for the two different (i.e., the up 
and down) spin directions [DODS, thus the designation 
"spin (s) unresricted"], and (ii) it relaxes the requirement 
that the electron orbitals be constrained by the symme- 
try of the external confining field [thus the designation 
"space (S) unrestricted"]. 

Subsequent to our discussion of symmetry breaking 
and formation of electron molecules in the framework 
of the sS-UHF method, we elaborate in section IV on 
the connection between the symmetry-broken sS-UHF 
solution and the exact one in the case of a clean 2c 
QD with B = 0. In particular, analysis of the exact 
solution using the conditional-probability-distribution 
(CPD) techniques reveals the formation of an electron 
molecule in agreement with the sS-UHF result. 

II. THOMSON QUANTUM DOT WITH OUTSIDE 
IMPURITIES 

To introduce some of the principal physical and 
methodological issues pertaining to symmetry breaking 
and formation of Wigner molecules in finite fermion sys- 
tems, we discuss first the 2e TQD at B = [we remind 
the reader that a sufficiently high value of Rw {= 2.39) 
was chosen, such that the QD is in the regime where the 
two electrons form a WM; for the other parameters used 
throughout this paper, see Ref. p3|] . The single-particle 
wave functions (modulus square) and total electron den- 
sities displayed in Fig. 1 are taken from calculations for 
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the 2e TQD in the presence of two attractive impurities 
(Qi = Qi = — le, represented in the figure by filled dots 
on the two sides of the QD) located symmetrically about 
the center of the QD at (x, y, z)= (±60, 0, 10) nm with 
the strength and location of the impurities purposefully 
chosen such that they will not affect the nature of the 
electronic ground states (for the same QD but without 
the impurities), except for orientational pinning in the 
case of symmetry-broken states (see below) . 

Constraining first the solution to maintain the sym- 
metry of the hamiltonian, in conjunction with double- 
occupancy of the HF orbitals by electrons of opposite 
spins [that is through the use of the RHF method with in- 
put trial electron densities satisfying the symmetry of the 
external potentials!!!!] , the resulting symmetry-adapted 
(SA) self-consistent singlet (<S) orbitals and correspond- 
ing total density distributions exhibit, as expected, an 
almost circular symmetry with minimal elliptical distor- 
tions (due to the impurities), see Fig. 1(a); without the 
impurities, the SA singlet is perfectly circularly symmet- 
ric. 

However, increasing the variational freedom through 
removal of the spatial symmetry and double-occupancy 
constraints via the use of the sS-UHF method results 
in a symmetry-broken singlet ground state of lower en- 
ergy, that is formation of a WM characterized by local- 
ized orbitals with the "bond length" (distance between 
the maxima in the total electronic distribution) equal to 
29 nm [see Fig. 1(b)]; the energy of this state is lower by 
1.62 meV than that of the SA solution [Fig. 1(a)]. This 
lowering of the ground-state energy reflects gain in cor- 
relation energy (for the definition of correlation energy 
see section IV below). Note that the 2e WM is orien- 
tationally pinned along the interimpurity axis. Similar 
formation of a WM occurs also for a 2e TQD without 
the pinning impurities (with an energy gain of 1.32 meV 
compared to the corresponding SA solution). However, 
in the absence of pinning, the formation_pf the WM is 
accompanied by orientational degeneracyE3 (that is there 
is an infinite manifold of rotationally degenerate sS-UHF 
ground states). 

The formation of a fcrmionic molecule, associated 
with electron localization, in the ground-state of a QD 
under magnetic-field-free conditions does not have an 
analog within the framework of the traditional models 
of atomic structure. u However, the physics underlying 
this phenomenon, which is a manifestation of SB re- 
sulting from strong electronic correlations (see also sec- 
tion IV below), resembles closely that found in doubly- 
excited two-electron atoms. Indeed, spectroscopical stud- 
ies on doubly-excited helium atoms revealed rovibra- 
tional bands which were interpreted, borrowing from 
models developed in the context of nuclear and molec- 
ular physics, by invoking the formation of a "triatomic" 
molecule comprised of the two localized electrons and the 
He 2+ nucleus (a particle), with the-iCollinear configura- 
tion being of particular significance. □ 

The first electronically excited state of the 2e TQD is 



the triplet T state (with the spins of the two electrons 
parallel to each other) whose total electron density dis- 
tribution [Fig. 1(c), left] resembles that of the ground- 
state singlet [Fig. 1(b), left]. However, the individual 
electronic wave functions in the T state differ in char- 
acter from those of the S state [compare right panels in 
Fig. 1(c) and 1(b)], with the lower-energy one being s- 
like (but elliptically distorted), and the other is a p-like 
orbital oriented by the impurities along the x-axis. Note 
that the T state has the symmetry of the hamiltonian 
including the two pinning impurities (that is, here the 
sS-UHF solution coincides with the symmetry-adapted 
one). The same orbital characters are obtained also in 
the absence of the pinning impurities, but without a pre- 
ferred orientation. In this case, however, the character of 
the T state is a result of spontaneous SB, with an energy 
gain of 0.09 meV compared to the corresponding SA (cir- 
cular) solution. Underlying the type of spontaneous SB 
in the (open-shell), T state of the 2e TQD is the Jahn- 
Teller (JT) effectEj where lowering of the total energy is 
achieved via mixing of the two- fold degenerate m = +1 
(p+) amd m — — 1 (p-) angular momentum states, con- 
comitant with a deformation of the_self-consistent poten- 
tial away from circular symmetryEH To distinguish such 
electron molecules from the WM discussed above for the 
closed-shell singlet state (whose formation is driven by 
the dominance of the electron-electron repulsion), we re- 
fer to them as JT electron molecules (JTEM's). Such 
spontaneous SB via JT distortions is familiar fromj-stud- 
ies of the rotational spectra of open-shell nucleiQa and 
from investigations of shape deformations of open-shell 
metal clusters olHI 

Similar calculations for the 2e TQD, but with repul- 
sive pinning impurities (that is Qi = Qi = +le) yield 
for the singlet ground state qualitatively similar results 
(with different values for the energies), but with an im- 
portant distinction that now the pinned orientation of the 
WM is rotated by tt/2 compared to the Q\ = Qi = — le 
case (i.e., the "intramolecular" axis of the WM is ori- 
ented normal to the interimpurity axis). Consequently, 
through variation of the sign (polarity) of the impurity 
gate voltages, one may "flip" the orientation of the WM, 
and with it the direction of the polarization of the elec- 
tronic charge distribution in the QD. In this way, the 
formation of WM's in QD's and the ability to control 
their discrete orientations via pinning voltage gates may 
serve as a method for the creation of <m-off information 
storage cells and nanoscale logic gates. EJ 



III. RUTHERFORD QUANTUM DOT 

Next, we examine the properties of a 2e Rutherford 
QD (RQD), that is a 2e QD with a central attractive im- 
purity. The sS-UHF singlet and triplet electronic orbitals 
corresponding to a 2e RQD for B = with a single impu- 
rity (Q = — 2e) located at (0,0,10) nm are shown in Fig. 
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2(b) and 2(c) respectively. They exhibit WM symme- 
try breaking and JT-distortion features similar to those 
found for the 2e TQD [compare Fig. l(b and c)], but with 
a reduced WM bond length and a more compact triplet. 
The "strength" of the SB depends of course on the impu- 
rity charge Q and/or its distance d from the QD plane. 
For example, for Q — — le and d = (and for an arbitary 
position of the impurity inside the QD), no symmetry 
breaking was found by us due to the strong trapping by 
the impurity of the two electrons which occupy circu- 
larly symmetric orbitals, resembling the behavior of the 
ground-state of the natural He atomic We also remark 
that for the case described in Fig. 2(b) and Fig. 2(c) both 
electrons are slightly trapped by the impurity potential 
for the S state, while for the T state the s-like electron is 
strongly trapped and the p-like electron occupies a much 
less bound orbital. 

Note that here, as with the case of a clean QD, the 
singlet WM and the JT-distorted triplet are free to rotate 
in the plane of the QD, since there are no off-centered 
pinning impurities. 

The large physical size of QD's makes them ideally 
suited for investigations of magnetic-field effects and con- 
trolled manipulations. To illustrate such effects, we dis- 
play in Fig. 2(a) the magnetic- field induced variation of 
the total energies of the S and T states in the 2e RQD 
(qualitatively similar behavior is found also for the 2e 
TQD). As seen, the energy of the singlet state increases 
and that of the triplet state decreases with increasing B; 
for fields B < B\ (0.2 T), the variation of the energy 
of the T (the slope of the curve) is smaller than that for 
B > B\ [see inset in Fig. 2(a)]. Furthermore, at a critical 
value B c = 2.8 T [marked by a down-arrow in Fig. 2(a)], 
the energies of the two states cross and from then on the 
triplet lies below the singlet. 

For the singlet, the broken-symmetry WM state [Fig. 
2(b)] maintains under the influence of the applied mag- 
netic field in the range considered in Fig. 2(a), with the 
increase in the magnetic-field strength leading to further 
shrinkage of the bond length of the WM accompanied by 
an overall increase of the energy of the S stateEj [see Fig. 
2(a)]. 

The influence of the magnetic field on the triplet state 
is more subtle. As aforementioned, at B = the sym- 
metry of the T state is broken by the JT effect involving 
mixing of the m = ±1 degenerate angular momentum 
orbitals (see Ref. [57|). On the other hand, the magnetic 
field lifts the degeneracy of these p + and p_ states (with- 
out mixing them), and this effect competes with the JT 
distortion. For small enough fields [B < Bi, see inset in 
Fig. 2(a)], the JT effect prevails, and thus the orbitals 
and electron densities remain similar to those shown in 
Fig. 2(c), and they maintain an orientational degeneracy 
in the plane. At stronger fields (B > Bi), the lifting of 
the energetic degeneracy of the p+ and p_ states over- 
comes the JT effect, and the second electron populates 
the lower of these two orbitals. As a result, circular sym- 
metry is recovered, as illustrated in Fig. 2(d). 



The overall decrease with B of the energy of the T 
state relative to the S state is due to enhanced stabiliza- 
tion by the (parallel-spin) exchange energy in the former, 
reduced Coulomb repulsion between the electrons occu- 
pying s- and p-like orbitals, and quenching of theJdnctic 
energy of the p-like orbital by the magnetic field.Ej This 
effect increases with B, and at B > B c the T state be- 
comes the ground state [see Fig. 2(a)]. Note that this 
transition is driven primarily by the interelectron repul- 
sion and not by the interaction of the electrons' mo- 
ments with the magnetic field (see Ref. || and references 
therein); for our system, the Zeeman splitting energy is 
0.026 mcV/T. At even larger fields, the Coulomb repul- 
sion between the electrons increases (due to the shrinking 
of the orbitals) resulting in an ascending trend of the en- 
ergy of the T state, which remains, however, lower than 
the singlet state. A similar scenario is found also for the 
TQD (without pinning impurities). 

In light of previous findingsEjO pertaining to forma- 
tion of fully spin-polarized symmetry-broken states in 
clean QD's (TQD's) at high magnetic fields, it is perti- 
nent to inquire whether the circular symmetry found for 
the T state of the 2e RQD for B 1 < B < 5 T will also be 
broken at higher fields. For the clean 2e TQD, we ver- 
ified that indeed an orientationally degenerate electron 
molecule [with the molecular orbitals of the electrons dis- 
tributed about the two molecular centers, see Fig. 2(e)] 
formed at sufficiently high B (e.g., B — 10 T). Such an 
electron molecule is akinc3 to the JTEM discussed above 
in the context of the triplet state for B = 0. Interest- 
ingly, such reemergence of a JTEM structure does not 
occur at these conditions for the 2e RQD studied here 
due to the enhanced gap between the p + orbital and the 
strongly trapped s orbital. This provides an additional 
venue for impurity-assisted manipulation and design of 
the electronic properties of QD's. 

IV. CONNECTION TO THE EXACT SOLUTION 

As mentioned earlier, for R\\r > 1, the sS-UHF ap- 
proach applied to QD's yields approximate ground-state 
solutions which violate the symmetries of the original 
many-body hamiltonian, e.g, the spontaneous breaking 
of rotational symmetry discussed in sections II and III 
for a circularly symmetric clean QD (i.e., a TQD) or one 
with a central impurity (i.e., RQD). At a first glance, 
this situation may appear puzzling, but it is not unique 
in the context of many-body theory of finite fcrmionic 
systems. Indeed, a similar situation was encountered 
in nuclear physics in the 1950's, when it was discov- 
ered that open shell nuclei carried permanent quadrupole 
moments and that many of them exhibited well devel- 
oped rotational spectra (i.e., they behaved like rigid ro- 
tors). The explanation of these experimental findings 
was formulated in the framework of breaking of the ro- 
tational symmetry associated with nuclear deformations 
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of the Jahn- Teller type, and it led to several celebrated 
models and semi-empirical methods, i.pr, the particle- 
plus-rotor model of Bohr and Mottelspa,cil the modified- 
anisotropic-oscillator modeljof NilssorO and the Strutin- 
sky shell-correction method.c3 In the language used by us 
here, this means that at the microscopic level the break- 
ing of the rotational symmetry had to be accounted^ for 
via space (S)-UHF methods (that is allowing the spatial- 
orbital solutions of the Hartree-Fock equations to assume 
symmetries lower than those of the underlying many- 
body hamiltonian). .— ■ 

Starting, with Peierls and Yoccozc3 (see also Peierls and 
Thoulcsalj), numerous theoretical investigationss have 
addressed the connection of the broken-symmetry HF 
solution to the exact solution, and led to the theora 
of restoration- of- symmetry via Projection Techniques .eZI 
The central physical ingredient of the Projection Tech- 
nique rests with the observation that a HF solution which 
breaks rotational symmetry is not unique, but belongs 
to an infinitely-degenerate manifold of states with dif- 
ferent spatial orientation. A proper linear combination 
of the HF determinants in such a manifold yields multi- 
determinental states with good total angular momenta 
that are a better approximation to the exact solution. 
For a comprehensive review on restoration of symmetry 
in the context of nuclear many-body theory, we refer the 
reader to the book by Ring and Schuck mentioned in Rcf. 
|35| (see also Ref. fl8], where the principles of restoration of 
rotational symmetry are discussed in the two-dimensional 
case). 

The restoration of symmetry via Projection Techniques 
in the case of sS-UHF solutions describing 2D-QD's will 
be presented in a future publication. In this section, we 
found it more convenient to discuss the connection be- 
tween the sS-UHF and the exact solutions by taking ad- 
vantage of the simplicity of solving the exact problem 
at B = for two electrons interacting via the Coulomb 
force and confined by an external parabolic confinement 
without impurities (clean QD). Indeed, it is well known 
that in this case, the exact Schrodinger equation for two 
interacting electrons is separable in the center-of-mass, 
R = (i"i + r 2 )/2 (with a corresponding mass M. = 2m*), 
and relative, r = ri — r 2 (with a corresponding reduced 
mass fi = m*/2), coordinates. 

To anal yse .the properties of the exact solutions, two 



quantities 



have customarily been extracted from the 



two-body wave function \I>(ri,r 2 ) [where the tilde de- 
notes both spatial and spin variables, i.e., r.^ = (r^s,) 
(i = 1, 2\r-see Appendix B]: (i) the usual pair-correlation 
functioned 



G(«)=27r(tf|*(ri-r 2 -v)|*) , 
and (ii) the electron density (ED) 

2 



l («) = (*|X)*(v-r i )|*) 



(3) 



(4) 



i=l 



However, for the exact "F in the case of a circularly 
symmetric confinement, both of these quantities turn out 
to be also circularly symmetric and thus they do not re- 
veal the full physical picture, the (possible) formation of 
an electron molecule generated by electron correlations. 
A more suitable quantity here is the conditional prob- 
ability distribution (CPD) "P(v|r 2 = Vo) for finding one 
electron at v given that the second electrpa is at r 2 = vo . 
This quantity has been extensively useaEI in the analy- 
sis of electron correlations in doubly-excited helium-like 
atoms, and is defined as follows, 

P(V ' r2 = V0) = <*|*(v -r 2 )|*> (5) 

[for details concerning the calculation of the quantities in 
Eqs. (|) - (|), see Appendix B]. 

In the left column of Fig. 3, we display the above three 
quantities for the exact ground state in the case of a 
parabolic confinement with the same parameters as used 
throughout this paper (see Ref. |23|) and B = 0. The 
pair-correlation function in Fig. 3(a) exhibits a well de- 
veloped depression at v = 0, namely the two electrons 
on the average keep apart from each other at a distance 
2r = 22.92 nm. In Fig. 3(b), we plot the CPD with 
v o = (fo, 0) (marked by a cross). It is seen that the max- 
imum probability for finding the second electron occurs 
at the diametrically opposite point — vo, in accordance 
with the picture of an electron molecule presented ear- 
lier within the sS-UHF approach. The exact electron 
density shown in Fig. 3(c) is circularly symmetric, as ex- 
pected. Comparing figures 3(b) and 3(c), the following 
interpretation ensues naturally, namely that the CPD in 
Fig. 3(b) describes the electron molecule in its intrinsic 
(body-fixed) frame of reference, while the electron den- 
sity in Fig. 3(c) describes the electron molecule in the 
laboratory frame of reference where rotational and and 
center-of-mass displacements are superimposed upon the 
intrinsic probability density. 

Fig. 3(d) displays the electron density for the corre- 
sponding sS-UHF ground state. As discussed earlier, the 
sS-UHF electron density breaks the rotational symme- 
try and clearly exhibits the morphology of an electron 
molecule, unlike the exact one in Fig. 3(c). It is apparent 
that the sS-UHF electron density corresponds to that in 
the intrinsic frame of the electron molecule. Restoration 
of the symmetry via Projection Techniques will bring the 
sS-UHF electron density closer to that of the exact solu- 
tion. As aforementioned, this interpretation is familiar in 
the context of nuclear and it is further supported by the 
CPD calculated with the sS-UHF ground-state [i.e., by 
using i]/ UHF in Eq. (||), instead of the exact many-body 
fy, see Appendix B] and displayed in Fig. 3(e). Although 
naturally not identical, the two CPD's (i.e., the exact 
and the UHF) are of similar nature and both illustrate 
graphically the correlation effect associated with electron 
localization and formation of an electron molecule. We 
further illustrate this point by contrasting the exact and 
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sS-UHF CPD's with the CPD of the symmetry-adapted 
RHF ground state, shown in Fig. 3(f). In this latter 
case, * RHF = ^ RHF (r 1 )^ RHF (r 2 ) x (si,s 2 - 1 S = 0), with 
■0 RHF (v) being the Is orbital of the RHF; as a result, 
the probability of finding one electron at v is independent 
of the position of the second electron and it is centered 
about the origin where it achieves its maximum value, as 
expected from an independent-particle description, i.e., 

P RHF (v|r a = v ) = |^ RHF (v)| 2 . (6) 

Finally it is of interest to examine the energetic as- 
pects of the symmetry breaking. Indeed the energy of 
the ground-state is E™ F = 22.77 meV for the RHF so- 
lution, £" HF = 21.45 meV for the sS-UHF solution, and 
E** = 19.80 meV for the exact one. Since the correlation 
energy is-the difference between the RHF and the exact 
energies]^ one finds e c g °" = £g HF ~ E™ = 2.97 meV. 
Thus the symmetry breaking associated with the sS-UHF 
solution is able to capture (E^ F - Eg HF )/e c g °" = u% q{ 
the correlation energy; the remaining amount can be cap=. 
tured through improvements via Projection techniques. til 

V. SUMMARY 

The dominance of interelectron interactions and cor- 
relations in quantum dots (often referred to as "arti- 
ficial atoms") results in spontaneous symmetry break- 
ing of the self-consistent mean field, beyond the central- 
field picture on which the modern understanding of 
atomic structure has been founded. Indeed, as shown in 
this Pfp^r through spin-and-space unrestricted Hartree- 
FocU13'l2I many-body calculations for the singlet ground 
state of a 2e QD, such symmetry breaking and the result- 
ing formation of Wigner molecules do occur for a wide 
range of system parameters (with or without impurities) 
even in the most elementary case of two-electron QD's. 

Additionally, in the case of fully polarized dots (e.g., 
the triplet state of a 2e QD), the sS-UHF calculations 
show that similar electron molecules form due to sym- 
metry breaking associated with Jahn- Teller-type distor- 
tions. 

In spite of the simplicity of the hamiltonian, correla- 
tions in two-electron QD's underlie a remarkably rich and 
complex physical behavior. As a result, 2e QD's can serve 
as paradigms for a unification of concepts pertaining to 
spontaneous SB in various finite fermion aggregates, in- 
cluding nuclei, excited atoms, and clusters. 

Furthermore, impurity and/or magnetic-field effects al- 
low for controlled manipulation of the formation and pin- 
ning of the discrete orientations of the electron molecules 
in 2e QD's. Such ability may be employed in future appli- 
cations of QD's as nanoscale logic cells and information 
storage elements. Impurity effects were also illustrated 
for the case of a quantum dot with more than two elec- 
trons. 



Confirmation of the formation in 2e QD's of electron 
molecules associated with symmetry breaking of the self- 
consistent field (in the context of sS-UHF calculations) 
was obtained through an analysis of the exact ground 
state via the conditional-probability-distribution tech- 
nique. 
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APPENDIX A: QD'S WITH MORE THAN TWO 
ELECTRONS 

In the main body of this paper we limited ourselves 
to the case of 2e QD's, since focusing on this elementary 
case allowed us to better elucidate the intricate physical 
principles involved in the formation and impurity control 
of electron molecules in QD's, from both the perspective 
of the sS-UHF treatment and the exact solution. In this 
appendix, we will present an example of the many differ- 
ent structural possibilities that can arise when impurities 
are introduced in a QD with a larger number of electrons. 
In particular, we consider a QD with 8 electrons in the 
presence of a hydrogenic-likc (d = 0) impurity of variable 
nominal charge Q — qe placed at the center (and in one 
case off-center) of the QD (for the other input parame- 
ters, we use same values as used throughout the text, see 
Ref. |23| ). Since this appendix does not intend to present 
an exhaustive study of larger dots, but simply aims at 
presenting an illustrative example, we will consider only 
one spin configuration, i.e., the sS-UHF solutions hav- 
ing 4 spin-up and 4 spin-down electrons. However, q will 
be allowed to take both positive (repulsive) and negative 
(attractive) values. 

Fig. 4(a) displays the sS-UHF electron density for 
a repulsive central impurity [located at (0,0,0)] with 
Q = +le. It is seen that a Wigner molecule consist- 
ing of a single ring with 8 electrons [denoted as a (0,8) 
ring] is formed in this case. For a slightly attractive cen- 
tral impurity with Q = — 0.2e, however, a structural 
change takes place, namely one electron moves to the 
center of the dot [see Fig. 4(b)]. Notice that this (1,7) 
structure agrees with the arrangement found in studies 
of classical noint-charges in a purely (Q — 0) harmonic 
confinement .Ej Increasing the attractive nominal charge 
to Q = — 0.5e [see Fig. 4(c)] results in a further structural 
change, namely the central impurity is now able to trap 
two electrons, thus leading to a (2,6) arrangement. A 
further increase of the attractive charge of the impurity 
to the value Q — — le does not produce any qualitative 
change in the (2,6) arrangement as long as the impurity 
remains at the center of the dot. By moving the impu- 
rity to an off-center position, however, various structural 
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morphologies can arise, an example of which is presented 
in Fig. 4(d) for a Q = — le impurity located at (20,0,0) 
nm, forming a highly distorted (2,6) WM (here the first 
index denotes that two electrons are trapped by the im- 
purity). 

A further increase in the attractive charge Q leads to 
sequential trapping of the remaining six electrons and to 
progressive elimination of symmetry breaking, until all 
eight electrons have been captured by the impurity (see 
also section III). 



APPENDIX B: TWO-ELECTRON WAVE 
FUNCTIONS AND THE DEFINITIONS (3)-(5) 

1. Exact solution for two electrons 

In the case of two interacting electrons confined by a 
parabolic potential (clean QD), one can perform a change 
of variables to center-of-mass (CM), R = (ri + r2)/2 and 
P = Pi + P2, and relative- motion (rm), r = ri r 2 and 
p = (pi — P2)/2, coordinates. Then the two-electron 
hamiltonian separates into CM and rm contributions, 



CM 



with 



and 



HcM = Hi + \ Muj2 ° r2 



p 2 1 1 o e 2 



2// 



2 2 

-uuj r 



(Bl) 



(B2) 



(B3) 



where M. = 2m*, R — |R|, fi = m*/2, and r = |r|. 

The center-of-mass motion associated with the coor- 
dinate R obeys a Schrodinger equation describing the 
motion of a particle of mass M. = 2m* in a 2D isotropic 
harmonic potential of frequency loq. Here ujq is the fre- 
quency of the original parabolic confinement, i.e., the in- 
terelectron repulsion has no baring on the center-of-mass 
motion. 

Using dimensionless polar coordinates U — R/(IqV2) 
and G, the center-of-mass wave function can be written 
as S(£/)e jMe with the radial part given by 

~(U) = C NM U\ M \e- u2 ^L^(UyCl) , (B4) 

where (N, M) are the radial and azimuthal (related 
to the angular momentum) quantum numbers, respec- 
tively, Co = 1/2, the normalization constant Cnm — 
[2iV!4l M l +1 /(7V + IMI)!] 1 / 2 , and l)™\x) are associated 
Laguerre polynomials. 

Since in the exact problem the Coulomb interaction 
preserves the rotational symmetry, the radial part of 
the wave function Q{u)e lm6 / \Ju associated with the 
relative motion obeys the following one-dimensional 



Schrodinger equation [in dimensionless polar coordinates 
u = r/{l V2) and 6], 



d 2 Q + ^-m 2 + 1/4 ^ RwV2 + £ 
du 2 u 2 u fru>o/2 



}n = . 
(B5) 



The 1/u term results from the interelectron repulsion. 

Defining U; = (u^s^) (where Si is the spin of the ith 
electron and i — 1,2), the exact many-body (here two- 
body) wave function is given by 

*(ui,u 2 ) = *(ui,u 2 )x(si,S2) 



1 n(u) 

27T y/u 



e* m "Z(Uy M v X (si,s 2 ) , (B6) 



with U; = rj/(Z v / 2) and x( s ii s 2) is the spin part. 
The exact eigenvalues are given by 

E NM ,nm = htoo{2N + \M\ + 1) + e(n, \m\) , (B7) 

where e(n, \m\) are the eigenvalues associated with the 
relative motion [see Eq. (pq)], (n, m) being the corre- 
sponding radial and azimuthal quantum numbers. 



2. Pair correlation, electron density, and conditional 
probability distribution for the exact and sS-UHF 
two-electron wave functions 

The bracket notation in Eqs. (||), (Q), and (||) imply 
integrations over both the spatial and spin variable s. I n 
the case of the exact wave function ^ given by Eq. flB(j ) , 
the spin variables separate out, and thus G(v), n(v), and 
the CPD can be expressed as double integrals over the 
positions (ri,r2) of the following two-body spatial prob- 
ability density VF(ri,r 2 ) = |$(ri, r 2 )| 2 , where $ is the 
spatial part of 'P. In particular, one finds 



G(v) =2tt J J 6(n - r 2 - v)W(r u r 2 )<fridr a . ( B 8) 

for the pair-correlation function, 

2 

n(v) = J J ^(v-rOWti-LraJdridra , (B9) 

i — 1 

for the electron density, and 

W(v,r 2 = v ) 



V{v\r 2 = v ) 



/ driW(ri,r 2 = v ) 



(BIO) 



for the conditional probability distribution. 

In the case of the sS-UHF singlet ground state, the two 
electrons occupy two different spatial orbitals -01 ( r ) an d 
tp2 (r). Then the corresponding two-body wave function 
is the following single determinant, 
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UHF 



(n,r 2 ) = -4[V'i(ri)a(l)V 2 (r 2 )/3(2) 



V2 



-Vi(r 2 )a(2)V 2 (ri)/3(l)] , (Bll) 



where a and /3 denote the spin- up and spin-down spinors, 
respectively. Integration of |^ UHF | 2 over the spin vari- 
ables yields the following two-body spatial probability 
density, 



W 



UHF 



(n,r 2 ) = -[\Mn)\ 2 \Mr2)\ 2 + |Vi(r 2 )| 2 |^(n)l 2 ] 

(B12) 



To calculate the Conditional Probability Distribution 
in the case of the sS-U HF ground state, one replaces W 
by W VHF in Eq. ( pTp| ). 



1 M.A. Kastner, Phys. Today 46(1), 24 (1993). 

2 R.C. Ashoori, Nature 379, 413 (1996). 

3 S. Tarucha, D.G. Austing, T. Honda, R.J. van der Hage, 
and L.P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 

4 E.U. Condon and G.H. Shortley, The theory of atomic spec- 
tra (Cambridge Univ. Press, London, 1935). 

5 D.R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). 

6 The failure of Bohr- type models, based on the orbiting 
of spatially correlated electrons in conjunction with the 
Bohr-Sommerfeld quantization rule, to yield a reasonable 
estimate of the ground state of the He atom signaled a 
looming crisis in physics in the 1920's, which Bohr himself, 
as well as others, had been keenly aware of, as summa- 
rized succintly by Sommerfeld: "All attempts made hith- 
erto to solve the problem of the neutral helium atom have 
proved to be unsuccessful" [A. Sommerfeld, Rev. Sci. Inst. 
7, 509 (1923)]; see also the 10th chapter entitled "It was 
the Spring of hope, it was the Winter of despair" in A. 
Pais, Niels-Bohr's Times (Clarendon, Oxford, 1991), the 
review by J.H. Van Vleck, Phil. Mag. 44, 842 (1922), 
and M. Born, The Mechanics of the Atom (Ungar, New 
York, 1927). While, since, numerical solutions of the two- 
electron Schrodinger equation provided a quantitative res- 
olution to the problem, the hrst successful semiclassical 
treatment of the three-body Coulomb system awaited till 
1980, see J.G. Leopold and I.C. Percival, J. Phys. B 13, 
1037 (1980); see also D. Wintgen, K. Richter and G. Tan- 
ner, Chaos 2, 19 (1992). Furthermore, based on rather gen- 
eral group-theoretical arguments arising from the observa- 
tion of hierarchies with lower symmetry in the excited spec- 
tra, and motivated by ideas originating in nuclear-physics 
spectroscopy, it has been discovered in the late 1970's and 
early 1980's that electron correlations in doubly excited 
He lead to quantization of the spectrum much like in a lin- 
ear triatomic molecule, e-He 2+ -e. This molecular picture, 
with near rigidity and separability, results in "infinite se- 
quences of vibrational levels, on each of which is built an in- 
finite sequence of rotational levels", see M.E. Kellman and 



D.R. Herrick, Phys. Rev. A 22, 1536 (1980); see also G.S. 
Ezra and R.S. Berry, Phys. Rev. A 28, 1974 (1983), and 
S. Watanabe and CD. Lin, Phys. Rev. A 34, 823 (1986). 
Interestingly, as in the aforementioned semiclassical treat- 
ments, the collinear configuration plays a special role in 
the molecule-like model, serving perhaps as "partial vindi- 
cation" of the geometry considered originally by Bohr. 

7 A. Bohr and B.R. Mottelson, Nuclear structure (Benjamin, 
Reading, MA, 1975). 

8 W. Nazarewicz, Nucl. Phys. A 574, 27c (1994). 

9 C. Yannouleas and U. Landman, Phys. Rev. B 51, 1902 

(1995) . 

10 W. Ekardt, Metal clusters Ch. 1 & Ch. 4 (Wiley, New York, 
1999). 

11 In modeling of practical "on-off" (i.e., two-state) devices, 
one will undoubtedly have to take into account additional 
potential terms in the many-body hamiltonian, such as 
those associated with image charges resulting from inter- 
faces between different dielectric materials. The feasibil- 
ity of constructing nanoscale logic cells via controlled ma- 
nipulation of the antipodal arrangements of two electrons 
in a larger device made out of four dots arranged in a 
square to form a quantum-dot supermolecule is currently 
under intense experimental investigation, see I. Amlani, 
A.O. Orlov, G. Toth, G.H. Bernstein, C.S. Lent, and G.L. 
Snider, Science 284, 289 (1999). 

12 E. Wigner, Phys. Rev. 46, 1002 (1934). 

13 P.A. Maksym, Phys. Rev. B 53, 10871 (1996). 

14 H.-M. Muller and S.E. Koonin, Phys. Rev. B 54, 14532 

(1996) . 

15 R. Egger, W. Hausler, CH. Mak, and H. Grabert, Phys. 
Rev. Lett. 82, 3320 (1999). 

16 C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 5325 
(1999). 

17 P.W. Anderson, Basic notions of condensed matter physics 
(Addison- Wesley, Reading, MA, 1984) . 

18 D.J. Thouless, The quantum mechanics of many-body sys- 
tems (Academic Press, New York, 1972). 

19 H. Fukutome, Int. J. Quant. Chem. 20, 955 (1981). 

20 For studies pertaining to the arrangement of classical point- 
charges in a harmonic confinement, see Yu.E. Lozovik, Usp. 
Fiz. Nauk 153, 356 (1987) [Sov. Phys. Usp. 30, 912 (1987)] 
and V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49, 
2667 (1994). 

21 That Rw is the only relevant parameter in the formation 
of a WM in a clean QD can be seen immediately from the 
exact Schrodinger equation describing the relative motion 
of two interacting electrons at B = in a parabolic confine- 
ment (Appendix B). Indeed, this equation when expressed 
in dimensionless variables as in Eq. (B5) contains Rw as 
the sole free parameter. The corresponding center-of-mass 
motion is independent of the interelectron repulsion. For 
how to extract information about the Wigner molecule 
from the exact two-body solution, see section IV. 

22 B. Tanatar and D.M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

23 In all calculations, we used p*(GaAs)= —0.44, m*(GaAs) = 
0.067m e , and huo = 5 meV, with Rw = 2.39 (that is k = 
8); the same value of Rw is obtained by using ft(GaAs) = 
12.9 and Tiujq = 1.9 meV, yielding the similar results given 



8 



here, but with different absoiute energy vaiues. 

24 M. Wagner, U. Merkt, and A.V. Chapiik, Phys. Rev. B 45, 
1951 (1992). 

25 D. Pfannkuche, V. Gudmundsson, and P.A. Maksym, Phys. 
Rev. B 47, 2244 (1993). 

26 M. Taut, J. Phys. A 27, 1045 (1994). 

27 M. Taut, Phys. Rev. A 48, 3561 (1993). 

28 M. Rontani, F. Rossi, F. Manghi, and E. Molinari, Phys. 
Rev. B 59, 10 165 (1999). 

29 Here we differentiate between the exact solution of the 
hamiltonian specified by Eqs. (|l|) and (^) from "exact diag- 
onalization methods" of Hubbard-type hamiltonians used 
as models of QD's [see e.g., R. Berkovits, Phys. Rev. Lett. 
81, 2128 (1998)]. 

30 The sS-UHF equations were solved using the Pople-Nesbet- 
Roothaan formalism, see p. 214 in Ref. [jlj; see also J. A. 
Pople and R.K. Nesbet, J. Chem. Phys. 22, 571 (1954). 
Apart from Ref. [l(], the only previous use of a Space-UHF 
method in QD's was in an investigation (see Ref. ^) of 
WM's but only for fully polarized QD's (i.e., under high 
magnetic fields where the spin unrestriction is not at play) . 

31 A. Szabo and N.S. Ostlund, Modern quantum chemistry 
(McGraw-Hill, New York, 1989). 

32 E.K.U. Gross, E. Runge, and O. Heinonen, Many-particle 
theory (Adam Hilger, Bristol, 1991) ch. 8. 

33 From the point of view of the terminology adopted here, 
the simple designation Hartree-Fock (HF) in the litera- 
ture most often refers to the RHF, in particular in atomic 
physics and the physics of the homogeneous electron gas. In 
nuclear physics, however, the simple designation HF most 
often refers to a space (S)-UHF. 

34 see the article by Ezra and Berry in Ref. [| see also R.S. 
Berry in Structure and dynamics of atoms and molecules: 
Conceptual trends, edited by J.L. Calais and E.S. Kryachko 
(Kluwer Academic, Dordrecht, 1995) p. 155, and references 
therein. 

35 The orientational degeneracy of the sS-UHF solutions de- 
scribing an electron molecule in a cirularly symmetric QD 
is a residue of the rotational invariance of the exact hamil- 
tonian and naturally it suggests the possible emergence of 
a ro-vibrational spectrum associated with the excitations 
of the electron molecule. Depending on the value of the 
parameter Rw, the electron molecule may exhibit differ- 
ent degrees of "rigidity", with electron molecules of high 
rigidity (forming for large Rw) exhibiting well-developed 
ro-vibrational bands; for the value of Rw = 2.39 used 
by us here, the electron mo lecule is rather "flo ppy" [C. 
Yannouleas and U. Landman, cond-mat/0003245|. Such ro- 
vibrational spectra can be described by several theoretical 
approaches, for example: (i) at a more phenomenological 
level, by a quantization of the collective classical motions 
associated with the molecule, such as rotations, and bend- 
ing and stretching vibratio ns [ see Kellman and Herrick in 
Ref. ^ and Maksym in Ref. [f3|]; (ii) at a microscopic level, 
by projecting the symmetry-broken HF determinant onto 
multi-determinental wave functions with good angular mo- 
menta, thus restoring the symmetry of the exact hamil- 
tonian; see P. Ring, and P. Schuck, The nuclear many- 
body problem, Ch. 11 ( Springer- Verlag, New York, 1980) 
and P.-O. Lowdin, in Quantum theory of atoms, molecules, 



and the solid state, Edited by P.-O. Lowdin, pp. 601 - 623 
(Academic Press, New York, 1966). Furthermore, starting 
from an exact solution for the electrons in the QD, rovi- 
brational characteristics can be obtained from the many- 
body eigenvalue spectra evaluated for different total angu- 
lar momenta and spins in conjunction with an analysis of 
the conditi onal- probability den sities [C. Yannouleas and U. 
Landman, cond-ma t/0003245| . We note that the discrete 
ro-vibrational collective spectra associated with symmetry- 
breaking in a QD may be viewed as finite analogs to the 
Goldstone modes accompanying symmetry breaking tran- 
sitions in extended media (see Ref. p^ ). 

36 H.A. Jahn and E. Teller, Proc. Roy. Soc. London Ser. A 
161, 220 (1937). 

37 Recall that for an electron moving in a 2D circularly sym- 
metric potential the wave functions are given (in polar coor- 
dinates) by ^n r ,m(r, 9) — R(r; n r , \m\) exp(im9) with cor- 
responding energies £(n r , \m\), where n r is the number of 
radial nodes and m is the angular momentum quantum 
number, with m = 0, ±1, ±2, ... In the JT framework, the 
second electron of the (open-shell) triplet state of the 2e 
QD occupies an orbital formed as a linear combination of 
the n r = 0, m = ±1 states (denoted by p+ and p~), i.e., 
p x cx p+ + p- or p y oc p+ — p_ , with the one of lower 
energy being selected through the self-consistent distortion 
of the central mean field. As a consequence, the orbitals of 
both electrons in the 2e QD exhibit noncircular symmetry. 
The energy gain due to the JT distortion depends on the 
strength of the confinement potential, and it increases for 
weaker confinements. In natural atoms, the spontaneous JT 
effect is fairly small and perturbative methods are usually 
adequate; see Ch. III.2 in Ref. 

38 This cancellation of crystallization due to a strong- 
attractive-impurity effect is illustrated in Fig. 4(d) where 
the attractive impurity is able to capture two of the eight 
electrons in the dot. The two captured electrons form a sin- 
glet subunit and their orbitals are clearly circularly sym- 
metric (see Appendix A below). 

39 The effective confinement frequency f2 e ff = (ujq + LJc/4) 1 ^ 2 
in the Darwin-Fock model of a 2D harmonically-confined 
electron in the presence of a magnetic field B (where wo is 
the frequency of the harmonic potential and u c = eB/m*c 
is the cyclotron frequency) increases as a function of the 
magnetic field. In a magnetic field, the energy of a harmon- 
ically confined electron is given in the Darwin-Fock model 
by En r , m = (2n r + \m\ + l)hQ eff - mhuj c /2. 

40 Such JT-like electron molecules have been called Wigner 
molecules in previous publications treating fully polarized 
QD's (Refs.|l|and|l|). However, this naming is not appro- 
priate since their formation is not driven by the dominance 
of the interelectron repulsion (over the zero-point energy), 
which for large ma gne tic fields can be treated perturba- 
tively (Refs. E4 and B6J) . 

41 A. Bohr and B.R. Mottelson, Mat. Fys. Medd. Dan. Vi- 
densk. Selsk. 27 (No. 16) (1953). 

42 S.G. Nilsson, Mat. Fys. Medd. Dan. Vidensk. Selsk. 29 
(No. 16) (1955). 

43 V.M. Strutinsky, Nucl. Phys. A95, 420 (1967). 

44 For an early review of the application of Hartree-Fock 
methods to the description of open-shell, deformed nuclei, 







see G. Ripka, Adv. Nucl. Phys. 1, 183 (1968). 

45 R.E. Peierls and J. Yoccoz, Proc. Phys. Soc. (London) 
A70, 381 (1957). 

46 R.E. Peierls and D.J. Thouless, Nucl. Phys. 38, 154 (1962). 

47 The concept of restoration of symmetry is not restricted 
to spatial symmetry, but includes also restoration of spin 
symmetry (see the paper by Lowdin in Ref. |35| 

48 G.E. Brown, Unified theory of models and forces (North- 
Holland, Amsterdam, 1967) p. 64. 

49 For the purely parabolic confinement considered r-kere, it 
turns out that the pair-correlation function reducesEj to the 
modulus square of the radial part of the wave function (see 
Appendix B) describing the relative motion, i.e., G(v) — 

n 2 (v)/v. 

50 See p. 620 of the article by Lowdin mentioned in Ref. ^ (in 
this article the RHF is referred to as conventional Hartree- 
Fock scheme); this definition was also followed in Ref. 
where a RHF solution was contrasted to the exact ground- 
state solution for two interacting electrons. 

51 C. Yannouleas and U. Landman, to be published. 

FIG. 1. Total electron densities (left frame in each pan- 
nel) and contours of modulus square of the individual orbitals 
(right frames) for the clean (Thomson) 2e TQD at B = 0, in 
the presence of two attractive (Qi = Q2 — — le) impurities 
whose projected positions on the QD plane (xy) are denoted 
by filled circles, (a) The symmetry-adapted (RHF) singlet 
(5). The energy of the SA singlet is higher by 1.34 meV than 
that of the corresponding triplet, manifesting a shortcoming 
of the RHF method, (b) The sS-UHF Wigner-molecule sin- 
glet, orientationally pinned along the inter-impurity axis and 
exhibiting localized orbitals, one to the left and the other to 
the right of the QD middle; note the lack of reflection sym- 
metry of the individual orbitals about the mirror (yz) plane 
normal to the inter-impurity axis and passing through the 
QD's center. When the WM singlet is rotated by n/2, the 
energy of the system rises by 0.58 meV. (c) The triplet (T) 
state with an elliptically deformed s-like and ap^-like orbital, 
orientationally pinned by the impurities. The energy of the 
symmetry-broken singlet in (b) is lower by 0.28 meV than 
that of the triplet in (c), compared to 0.43 meV for the TQD 
without pinning impurities. Distances are in nm and the elec- 
tron densities in 10 -4 nm -2 . The parameters characterizing 
the QD shown here are those given in Ref. E3l 



FIG. 2. (a-d) Energetics and individual orbitals (modu- 
lus square) from sS-UHF calculations for the 2e (Rutherford) 
RQD, with a central impurity Q — — 2e at (0,0,10) nm. En- 
ergies (in meV) of the S and T states versus B (in Tesla), 
crossing at B c = 2.8 T (marked by an arrow), are shown 
in (a). An expanded view of the energy of the T state for 
small fields, exhibiting a transition from the spontaneous JT 
regime to a circular symmetric state at B > B\ (B\= 0.2 T), 
is shown in the inset; included also is the energy of the SA 
solution (dashed curve). At B = 0, contours of the orbitals 
of the two electrons for the WM singlet are shown (superim- 
posed on each other) in (b), with one of the orbitals depicted 
by a solid line and the other by a dashed line. The orbitals for 
the spontaneously JT-distorted triplet at B = are shown in 
the two panels of (c). For the singlet, the energy gain due to 
SB (that is lowering of the total energy with respect to the 
symmetry-adapted RHF state) is 0.39 meV and the energy 
gain due to the JT-distortion of the triplet is 0.07 meV; the 
energy difference between the T and S states is 2.01 meV. 
The transition to a circular symmetric triplet is illustrated 
for B = 0.6 T in (d). (e) Individual orbitals for the sym- 
metry-broken ground-state triplet of the clean (Thomson) 2e 
TQD at B = 10 T. Under the same conditions, the ground 
state triplet in the 2e RQD remains circularly symmetric. Dis- 
tances in nm and orbital densities in 10 -4 nm~ 2 . Note the 
different length scales of the (x, y) axes in (b-e) compared 
to those in Fig. 1. The parameters characterizing the QD's 
shown here are those given in Ref. ^3| 

FIG. 3. Ground-state results for a clean 2e QD with the 
parameters given in Ref. ^ and B — 0; in all cases the ground 
state is a singlet, (a-c) correspond to the exact solution, (a) 
The pair correlation function [G(v), see Eq. (H)] plotted ver- 
sus v exhibiting a well developed depression at v = 0, with 
a mean electron-to-electron separation of 2ro = 22.92 nm. 
(b) The electron conditional probability distribution [CPD, 
see Eq. (Q)] with vo = (ro,Q) (denoted by a cross), showing 
formation of a 2e electron molecule, (c) The electron den- 
sity [ED, see Eq. (^J)], reflecting the conservation of circular 
symmetry by the exact two-electron solution. The results 
displayed in (d) and (e) correspond to calculations using the 
sS-UHF method. Spontaneous symmetry breaking leading to 
formation of a 2e-molecule (with a bond length of 2r' = 28.16 
nm) is exhibited by the sS-UHF electron density shown in 
(d). (e) Such an electron molecule is further reflected in the 
sS-UHF CPD with v = (r' , 0) (denoted by a cross), (f) The 
CPD with vo = (r'o,0) (denoted by a cross) corresponding 
to the restricted Hartree-Fock (RHF) ground state (i.e., the 
symmetry-adapted solution) exhibits, as expected, no sym- 
metry-breaking signature. Lengths are in nm and density 
functions [G(v), ED's, and CPD's] are in units of 10 -4 nm -2 . 
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FIG. 4. Electron densities obtained via sS-UHF calcula- 
tions of a QD containing 8 electrons at B = in the presence 
of a central [located at (0,0,0) in (a-c)] and an off-centerd 
[located at (20,0,0) nm in (d)] hydrogenic impurity. The pa- 
rameters characterizing the QD shown here are those given 
in Ref. ^3] and the charge Q of the impurity is as marked 
on the figure. All the cases correspond to zero total spin, 
i.e., 4e with spin up and 4e with spin down. For each of the 
cases, we observe formation of a Wigner molecule, with its 
structure dependent on the polarity (sign), magnitude and 
location of the impuriry. (a) A repulsive central impurity 
(Q = +le), resulting in an 8e WM with a ring structure 
and an empty-electron central region, denoted as (0,8). (b) 
Slightly attractive central impurity (Q = — 0.2e) leading to 
formation of an (1,7) two-ring WM. (c) A stronger attractive 
central impurity (Q — — 0.5e) associated with a (2,6) WM. 
(d) An off-center attractive impurity (Q = — le), resulting in 
a distorted 6e WM with two electrons trapped by the impu- 
rity. Lengths in nm and electron densities in units of 10 -3 
nm -2 . 
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